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We introduce a family of variational ansatz states for chains of anyons which optimally exploits 
the structure of the anyonic Hilbert space. This ansatz is the natural analog of the multi-scale 
entanglement renormalization ansatz for spin chains. In particular, it has the same interpretation as 
a coarse-graining procedure and is expected to accurately describe critical systems with algebraically 
decaying correlations. We numerically investigate the validity of this ansatz using the anyonic golden 
chain and its relatives as a testbed. This demonstrates the power of entanglement renormalization 
in a setting with non-abelian exchange statistics, extending previous work on qudits, bosons and 
fermions in two dimensions. 



I. INTRODUCTION 

Strongly correlated quantum many-body systems with 
topological order have been proposed as a substrate 
for building fault-tolerant quantum computers Un- 
der this proposal, logical information is stored within a 
subspace of a fixed number of quasi-particles (anyons). 
Computation is performed by exploiting the non-abelian 
statistics obeyed by these anyons under braiding. Com- 
pared to more conventional implementations of quantum 
computers, this offers an intrinsic resilience to noise: lo- 
cal perturbations cannot decohere the stored information 
because of the non-local nature of the encoding. 

At present, perhaps the most promising candidate 
systems exhibiting non-abelian statistics are fractional 
quantum Hall systems^"''. Some amount of experimental 
evidence is already available in this setting^"^^. Other 
proposed systems include topological insulators^^"^* and 
lattice spin systems^'^^'^^ whose local interactions could 
be enginereed artificially^^'^®. Independently of the pro- 
posed physical realization, the stability of the topologi- 
cally ordered phase with respect to (local) perturbations 
is of great interest for topological quantum computation. 
This presents a formidable theoretical challenge. Possi- 
ble approaches range from the study of concrete physical 
models (see e.g.,^^"^^) or bounds on the gap for general 
families of models^'^"^^ to the investigation of effective 
Hamiltonians describing the inter-anyon interactions (see 
e.g.,^^). In the context of the latter, it is natural to 
study paradigmatic models such as 1-dimensional any- 
onic chains. These are the natural counterpart of spin 
chains. The study of such systems has led, e.g., to an ex- 
actly solved model called the golden chain^^, and novel 
realizations of infinite randomness critical phases^^"'^^. 

A peculiar aspect of anyonic systems is the structure of 
the Hilbert space of N anyons. In contrast to the space 
of N qudits, this space does not decompose into an A^- 
fold tensor product. Instead, its dimension scales as , 
where the quantum dimension d is generally not an in- 
teger. For analytical studies, it is sometimes convenient 
to embed this space into a larger tensor product space. 
However, this approach is inconvenient when using varia- 
tional methods: with the most straightforward encoding, 
it may be unclear how to vary (numerically) over the 



subset of physical states. Furthermore, this leads to a 
significant increase in computational complexity, particu- 
larly because locality is only approximately preserved by 
such an embedding. Even for one-dimensional anyonic 
systems, these issues complicate the direct application 
numerical methods such as DMRG'^-'^. While such meth- 
ods may presently be the most powerful and successful 
tools available, these difficulties motivate the search for 
alternative approaches. 

There is another more profound reason why the naive 
application of numerical methods for qudits may be sub- 
optimal: local anyonic operators preserve the total topo- 
logical charge in their support. Since a realisitic Hamil- 
tonian consists of locally acting terms, topological charge 
conservation severly constrains the action of the Hamil- 
tonian both at the local level as well as on larger scales. 
Incorporating this fact into the method of study should 
therefore be highly beneficial. With a qudit embedding, 
the meaning of this constraint is obscured and may be 
hard to make use of. 

The purpose of this paper is to introduce a variational 
method for anyonic systems which avoids using unphys- 
ical additional states and optimally exploits the special 
structure of the Hilbert space. This is motivated in part 
by the goal of facilitating numerical studies. Perhaps 
more importantly, variational families of ansatz states 
can provide significant insight into the nature of quan- 
tum correlations in a given system. Even in this regard, 
a qudit embedding of anyons is generally undesirable, as 
a physically motivated ansatz for qudit systems may lose 
its significance when applied directly to the anyonic set- 
ting. 

Our scheme is inspired by renormalization group stud- 
ies of anyonic chains^*"'^° , composite anyon distillation'^^ 
and the multi-scale entanglement renormalization ansatz 
(MERA)'''^''^'' for (non-anyonic) spin chains. In fact, it 
can be understood as the natural anyonic counterpart of 
the latter and shares many of its properties. In partic- 
ular, it can be seen as a renormalization group scheme 
and is thus especially suited for describing scale-invariant 
systems. The scheme allows - in principle - to extract, 
e.g., critical exponents when considering fixed points of 
the renormalization group flow. More generally, it pro- 
vides procedures for computing parameters of the corre- 
sponding conformal field theory (CFT) in the continuum 
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limit. While such procedures are also available for other 
variational methods (e.g., matrix product states'^'"'"'^^ . by 
transfer matrix methods), they are particularly natural 
in the present context due to the scale-invariant form of 
the ansatz. 

Wc formulate our ansatz for one-dimensional (peri- 
odic) chains of anyons. This allows us to numerically 
test its validity for the golden chain^^ and its relatives^". 
Ultimately, though, it is desirable to find methods for 
two-dimensional systems. Our work makes some progress 
towards this goal: the non-anyonic MERA extends nat- 
urally to two dimensions and this is also the case for its 
anyonic counterpart. In the appendix, we bricfiy com- 
ment on such generalizations. 

Entanglement renormalization, while originally defined 
for qudit systems, has been extended to free bosons^^ and 
interacting fermions*^. Anyonic statistics are peculiar to 
two dimensions, encompassing both fermions and bosons 
(which are abelian) in addition to interesting non-abelian 
generalizations. Our ansatz applies to all such models, 
but its use is especially suggestive in the non-abelian case. 
This is due to the special form of the Hilbert space men- 
tioned earlier. Interestingly, the exact details of the ex- 
change statistics play essentially no role in the definition 
of the ansatz, though they become important when eval- 
uating expectation values of local operators in the case 
of e.g., two-dimensional arrangements of anyons. 

In a wider context, entanglement renormalization is 
a special instance of the class of tensor network states, 
which also includes e.g., MPS^^"^^, PEPS^^ as weU as the 
closely related TERG-states^^. These have wide appli- 
cability beyond qudit systems. In particular, a general 
framework for fermions has been developed^^^^^ . Our 
focus is on entanglement renormalization because of its 
unique operational interpretation, as well as the possi- 
bility of efficiently computing expectation values with- 
out approximations, even in two dimensions (in contrast 
to e.g., PEPS^°). Furthermore, anyonic entanglement 
renormalization is conceptually related to previous ana- 
lytical studies for anyonic chains^*"'^^. However, our work 
suggests that anyonic generalizations of other tensor net- 
work states should also be possible along similar lines, 
though further work may be needed to evaluate their de- 
scriptive power. 

The structure of this paper is as follows. In Section II, 
we review the formulation of the MERA for spin chains. 
In Section III, we give some background on anyonic sys- 
tems and their description in terms of fusion diagrams. 
We then present the anyonic entanglement renormaliza- 
tion ansatz in Section IV and show how to efficiently 
evaluate expectation values of local observables and cor- 
relation functions. We also discuss an operational inter- 
pretation in the context of composite anyon distillation. 
In Section V, we apply the ansatz to the golden chain 
and identify a renormalization group fixed point. Wc 
conclude in Section VI. 



II. MULTI-SCALE ENTANGLEMENT 
RENORMALIZATION 

To motivate our ansatz, we recapitulate the definition 
of the MERA for ID qudit chains and that of related ten- 
sor network states'^'^^'^'*. For more details, we refer to'"'^. 
The assumption underlying this ansatz is that certain lo- 
cal degrees of freedom are insignificant for the long-range 
character of the target state. The MERA is a quantum 
circuit which completely decouples and discards these de- 
grees of freedom using local unitaries. This procedure 
is subsequently repeated for the resulting coarse-grained 
description of the state. As a variational family of states, 
the MERA consists of those states which are turned into 
(A'^-fold) product states by such quantum circuits. 

It is convenient to use the diagrammatic formalism of 
tensor networks to describe the unitaries and the local 
isometrics constituting a MERA (the latter correspond 
to the elimination of local degrees of freedom, but the 
division into "disentangling" unitaries and isometrics is 
often arbitrary) . In this formalism, the identity on C' is 
represented by a single directed edge, and operators O : 
(•£.rf-j«)n _^ (^£d^(g)m represented by shaded boxes with 
n ingoing and m outgoing (ordered) edges. Labels on the 
edges from 1, . . . , d correspond to the elements of a fixed 
orthonormal basis of C^. Products of operators are taken 
by connecting outgoing with ingoing edges, and summing 
over labels of edges with no free ends. The trace of an 
operator on (C*)**" is taken by connecting each outgoing 
strand with the corresponding ingoing strand, and then 
contracting the tensor network (i.e., summing over all 
edge labelings) . Partial traces are computed analogously 
by connecting up subsets of edges. Tensor products are 
obtained by placing diagrams next to each other. 



With these conventions, the property W^^W = I 



of an isometry W : (C^)^" 
the simple form 



(C'i)« 



(C^)®™ (for m > n) takes 
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Identity (1) is crucial for the definition of the MERA 
ansatz. 

Consider a system of qudits arranged on a line with 
periodic boundary condtions. For any tensor network 
as shown in Fig. 1, one obtains a variational family of 
states l^f) = |^'|tp),{w}) G (C^)*^^ parametrized by the 
isometrics {W} and the state \(p) E at the top of 
the structure. This is the MERA ansatz. Fig. 1 repre- 
sents a map (C'')'^^ — >■ C which can be understood as 
the bra of the represented state. The fact that the 
recipe specifies instead of \'^) is a matter of prefer- 
ence, but the chosen convention has a natural operational 
interpretation: One may think of Fig. 1 as a renormal- 
ization prescription by decomposing the tensor network 
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(a) 



FIG. 1: Two examples of MERA structures: For a periodic 
qudit chain (at the bottom of the figure), a variational family 
of states is obtained by varying over the contents of the boxes. 
In (a), these are (adjoints of) isometries W : (C*)®^ C*, 
unitaries U : (C'')®^ _^ (C'')®^ and (the adjoint of) a state 
\(p) G C* for the strip at the top. Tree-like structures may also 
be considered, and lead to particularly simple expressions in 
the homogenous case'''^. Note that we choose to represent 
the top tensor by a (periodic) strip; the reason for this will 
become obvious once we move to the anyonic setting. 



along difTerent horizontal cuts Li. The strips [Li,Li+i[ 
represent coarse-graining maps which reduce the number 
of degrees of freedom at each level. 

Importantly, local expectation values and pair corre- 
lation functions can be computed efficiently for such a 
state. This is because the expectation value of an observ- 
able O is given by the contraction of the tensor network 
obtained by sandwiching O between the MERA network 
and its adjoint as in Fig. 2 (b). Eq. (1) then allows to 
simplify the network resulting in a significantly smaller 
network corresponding to the "causal cone" of the op- 
erator, see Fig. 2 (c). This network can be efficiently 
contracted, as the number of tensors scales logarithmi- 
cally with the number of sites N. For example, for the 
case of Fig. 2 (b), we get 



{^0\^) = o . . . £l„ o Tl„{OM (2) 

where Lmax = log2 — 1 and the superopera- 
tors associated to a level L are 
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FIG. 2: Tensor networks associated with the MERA of 
Fig. 1 (a). The network (a) represents the projection |^')(^' 
onto the state described by the MERA. The contraction of 
the network (b) gives the expectation value (4'|0|4') of a lo- 
cal operator O. Fig. (c) gives the same value as (b), and is 
obtained from it by using (1) repeatedly. 



for L < Lmax and 



O 



VfT/ 
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o 
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(4) 



In particular, the reduced two-site density operator of 
the state j^*) can be computed using the adjoint super- 
operators as 



tr„_2 I*) (1-1 = J^l o£l-- o £i_(|v.)(^|) . 

Similar expressions can be obtained for translates as well 
as for pair correlation functions evaluated at certain spe- 
cific distances^^"^"*. 

Given a MERA-structure as in Fig. 1, we have de- 
scribed a recipe giving a variational family of states for a 
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chain of N d-dimcnsional qudits parametrized by isome- 
tries of the form W : (C*)®" ^ (^d)®™ (the con- 
stants n > m depend on the MERA-structure). There 
is a natural way of enlarging this family: the edges (in 
higher levels of the network) may be interpreted as corre- 
sponding to C-^, where the refinement parameter x > d 
is larger than the dimension d of the physical qudits. 
Increasing this so-called bond dimension x amounts to 
using isometrics W : (C-*^)®" (cx)«"" (and corre- 
spondingly a state \(p) e C-^ at the top). To motivate 
a similar refinement in the anyonic setting, we point out 
that if X = is an integer power of the physical qu- 
dit dimension d, then this enlargement of the variational 
family is equivalent to replacing each edge in the MERA- 
structure by s edges. Note also that, since the number 
of isometrics is of order O(A^logA^) and each isometry 
is described by fewer than ^'^+™ parameters, the total 
number of parameters describing a MERA is of order 
Oipoly{x)NlogN). 

We conclude this short exposition by mentioning two 
subclasses of MERA-ansatz states which are of partic- 
ular interest: a translation-invariant MERA has iden- 
tical isometries within every level. As a consequence, 
the number of parameters of such a MERA scales 
as O {poly (x) log N). A scale-invariant MERA is one 
where all isometrics (with identical domain and range) 
arc chosen to be the same: Here the number of param- 
eters is only 0{poly{x)) independent of (see^"* for a 
more detailed discussion of the complexity of MERA 
and^^ for concrete examples). For a scale-invariant 
MERA, it is possible to write down two-point and three- 
point correlators explicitly in terms of eigenvalues of a 
certain coarse-graining superoperator defined in terms of 
the isometries as in (3). This allows to numerically ex- 
tract critical exponents and paramters from the associ- 
ated GET, as described e.g., in^^"^"*. 



their physical realization, their state space and exchange 
statistics are described by the axioms of a topological 
quantum field theory (TQFT) (see e.g.,^^). This formal- 
ism is extraordinarily useful for studying low-energy pro- 
cesses^^ as well as for the application to quantum compu- 
tation^, as it abstracts out the relevant physics: it spec- 
ifics the particle content, i.e., what particle types occur, 
describes what their internal degrees of freedom are, and 
how they are affected by braiding (exchanges) and fusion 
(which corresponds to bringing particles together). The 
relation between states of (and operations on) the phys- 
ical system and the abstract anyonic state space is dis- 
cussed extensively in the literature, see e.g.,^'^'^'^^'^^'^^. 

The algebraic object underlying such an anyonic theory 
is a modular category. Roughly, this consists of (i) a finite 
set of particle types f2 equipped with an involution * : 
51 — > f2 and containing a distinguished trivial particle 1 € 
57, (ii) fusion rules, i.e., a set of allowed triples of particles, 
(iii) a quantum dimension da > associated to every 
particle a, (iv) a tensor F indexed by 6 particles and (v) a 
3-index tensor R. These are required to satisfy a number 
of consistency conditions (see e.g.,^) the most important 
of which express associativity of fusion and compatibility 
of fusion with braiding. 

In more physical terms, the involution associates an 
antiparticle a* to every particle a, with 1* = 1 corre- 
sponding to the absence of a particle. The fusion rule 
summarizes the possible outcomes when bringing two 
particles together. The quantum dimensions give a rough 
measure of the growth of the anyonic state space when 
adding particles, and the F-tensor relates different bases 
of this space (as explained below). Finally, the i?-tensor 
encodes braiding of pairs of particles. This will be mostly 
irrelevant for our discussion, but becomes relevant for 2- 
dimensional arrangements of anyons as explained in Ap- 
pendix B. 



III. ANYONIC STATES AND OPERATORS 

In this section, we provide a short introduction to 
anyons, emphasizing the aspects relevant to the defini- 
tion of the anyonic entanglement renormalization ansatz: 
the anyonic Hilbert space and the isotopy-invariant for- 
malism for describing states and operators. We also dis- 
cuss the origin of the dynamics, that is, the definition of 
Hamiltonians for anyons. For a thorough and accessible 
introduction to anyons and topological quantum compu- 
tation, we recommend^ (see also^ for a recent review). 

A. A unified treatment of topological order 

Anyons arise as localized quasi-particle excitations in 
what can roughly be referred to as two-dimensional topo- 
logically ordered quantum media, e.g., qudit lattice sys- 
tems with certain Hamiltonians^'^^'^^, quantum Hall sys- 
tems^"^ or topological insulators^^"^^. Independently of 



B. The anyonic Hilbert space and anyon diagrams 

We proceed by explaining the construction of the any- 
onic Hilbert space using the data (i)-(iv) of a modular 
category. We will mostly follow the detailed exposition 
in^^, but will require slightly more general definitions 
when dealing with operators (see e.g., appendix of^^ for 
more details). 

The state space of a set of anyons depends on their 
types and on the surface the quantum medium is em- 
bedded in. We will discuss two cases in detail below: 
anyons pinned to fixed locations on a disc and on a 
torus. Usually, we assume that these are arranged on 
a chain, though one may also consider e.g., regular two- 
dimensional lattices; their particular geometric arrange- 
ment only becomes important when considering Hamil- 
tonians, but does not affect the definition of the Hilbert 
space. 

Starting point are certain trivalent graphs with di- 
rected edges. These correspond to pants decomposition 
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of the surface with punctures inserted at the anyons' posi- 
tions. Fixing such a graph, a basis of the anyonic Hilbert 
space is given by ah labehngs of the edges with particle 
labels from n satisfying the fusion rules at every vertex. 
We will give explicit examples for the punctured sphere 
and the torus below. 

Labeled graphs related by reversing the direction of 
an edge and simultaneously replacing its particle label 
by the associated antiparticle represent the same vectors. 
This is somewhat analogous to the formalism of Feynman 
diagrams. Indeed, anyon diagrams may, to some extent, 
be interpreted as particle world-lines, though this anal- 
ogy has its limitations. Note also that in many anyon the- 
ories of interest, such as the Fibonacci category consid- 
ered below, every particle is its own antiparticle, a* = a, 
and it is sufficient to work with undirected graphs. 

Dividing up trivalent graphs into neighborhoods of 
their vertices, one arrives at the following alternative de- 
scription: the total anyonic Hilbert space is the direct 
sum of tensor products of two-anyon fusion spaces V^^^ 
(respectively their dual splitting space V^'') correspond- 
ing to every vertex, where a,h,c € f2, and where the 
sum is taken over all fusion-consistent labelings (see be- 
low). The space can be thought of as the internal 
degrees of freedom of two anyons of type a and h whose 
combined topological charge is c. Equivalently, it is the 
space of two anyons a and 5 on a disc with total topo- 
logical charge c at the boundary. The latter is - in prin- 
ciple - a measurable quantity. We assume for simplicity 
that fusion is multiplicity-free, i.e., dimV^'^j, G {0, 1}, but 
our techniques directly generalize to models with fusion 
multiplicities (see^*" for the necessary adaptations in the 
diagrammatic calculus) . 

We pick a normalized vector \ah\ c) £ V^^ in each split- 
ting space and represent it using the isotopy-invariant 
formalism^^ as 



\ah\ 



{ab; c\ ~ 



dadb 



dadh 



1/4 



1/4 



ab 




where a scalar da called quantum dimension is associated 
to every label a. More generally, going from a vector to 
its dual vector corresponds to flipping the diagram along 
the horizontal axis and reversing all the arrows. Isotopy 
invariance means that diagrams may be continuously de- 
formed as long as endpoints are held fixed and edges are 
not passed through each other or around open endpoints. 



standard basis given by the vectors 

cti a2 as ffln 



dc 



1/4 




(5) 



where 6 G fi" ^ is such that this diagram is fusion- 
consistent. This corresponds to the decomposition 



One may switch between different bases using the (uni- 
tary) F move, i.e., the isomorphism between the two de- 
compositions of 

yabc ^ yab ^ yec ^ ybc ^ yaf 

f 

which is specified by the coefficients [-F'j''"^]e,/ in 
a b c a b c 



T Jd 



(6) 



The matrices [F^^'^] collectively constitute the i^-tensor 
of the modular category mentioned earlier. Matrix ele- 
ments of related isomorphisms such as 

= ® - v/, ® vf 
f 

can be expressed as functions of these. Diagrammatically, 
all basis changes can be computed using (6), the identity 



a{\b 



dndh 



(7) 



and the convention that lines with the trivial label may 
be added and removed arbitrarily, i.e.. 



(8) 




1. Anyons on a disc 

Consider the space of n anyons of types a G fi" 
pinned to fixed locations on a disc with total charge c G 51 
at the boundary. For this space, we will often use the 



An operator US' (c) : — taking the fusion space 
of TO anyons of types a = (ai, . . . , am) G fl™ with total 
charge c G to n anyons of types a' — {a'^, . . . , aJJ G il" 
with total charge c can be represented as a linear com- 
bination of trivalent graphs with to ingoing edges and 
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n outgoing edges attached to open endpoints and car- 
rying the corresponding labels. For example, using the 
standard basis (5), such an operators can be written as 



Ui (c) = T[US' {c)]g -,\a', 9; c)nm{S, b; c\ 



b,b' 



(9) 




b,n-2 



(10) 



ffli a2 as 
with the normalization factor 

1 




Any operator US' : V^f — > 0^, V^' taking m anyons 
of types a = (oi, . . . , a,„) £ fi™ to n anyons of type a' = 
{a[, . . . ,aJJ e fi" can also be represented in this fash- 
ion if it conserves the total charge, i.e., if it has block- 
diagonal form 



(11) 



Such operators act locally on subsets of n anyons of spec- 
ified types a € f2" and map them to a state of anyons of 
specified types a' € fi™. In fact, (10) (respectively (11)) 
represent the most general locally acting operator, their 
single most important property being charge preserva- 
tion. Additional properties such as unitarity impose fur- 
ther conditions on the form of the map / matrix US (c) for 
each c. 

More generally, we are interested in operators which 
act between spaces of the form 0^ (05:(Er Kf)i where 
r„ C il" specifies a set of n-tuples of available (spa- 
tial) particle configurations on a chain. A total charge- 
preserving operator of this kind takes the form 



(12) 



where t/^:(c) e End(e3^p^K",ea'6r„K°')- ™' 
can be understood as the projection of a general charge- 
preserving operator onto inputs and outputs from Tm 
and r„, respectively. More precisely, let 



denote the projection onto the subset of states with anyon 
labels from r„ defined by 



^ E \a,b;c),in{a,b;c\ 

E 



ben"" 



3er„ 



02 



as 



Then 



(13) 



which expresses the fact that the domain and range of 
the operator is restricted to states with certain particle 
configurations. 

Similarly as for tensor network states, we represent 
charge-preserving operators of the form (12) taking the 
space of m anyons to the space of n anyons by shaded 
boxes with m ingoing and n outgoing (unlabeled) edges; 
summation over these edge labels is left implicit. It is 
important to note, however, that such a box represents a 
different object compared to the case of tensor networks: 
it is defined by a family of maps {C^p™ (c)}c(eO which spec- 
ify a weighted superposition of certain labeled trivalent 
graphs embedded in the box, with m in- and n outgoing 
edges as in (10). 

The diagrammatic representation of anyonic operators 
and vectors satisfies simple rules with respect to com- 
position, tensor products and (partial) tracing. Adjoint 
operators arc obtained by flipping the diagram, reversing 
the arrows and complex conjugating all coefficients. Op- 
erators are multiplied or applied to vectors by stacking 
their representations on top of each other and connecting 
up out- with ingoing edges. Tensor products are obtained 
by placing diagrams next to each other. Traces and par- 
tial traces are the result of connecting up in- and out- 
going strands of an operator, with an additional minor 
modification (see Appendix A 2). 

These rules are remarkably similar to the contraction 
of tensor networks, although their origin is distinct. In 
particular, the notion of evaluating of a diagram is very 
different: A tensor network associates a scalar quantity 
to every labeling of the edges. This means that contrac- 
tion, i.e., summing over all labelings, results in a scalar. 
In contrast, anyonic diagrams associate a trivalent graph 
to every labeling, and the contraction results in a for- 
mal linear combination of equivalence classes of labeled 
trivalent graphs. Equivalence is defined by isotopy and 
the local rules (6), (7) and (8). 



2. Anyons on a torus 

Here we are interested in periodic chains of anyons ar- 
ranged on a line, and a few modifications of the above 
formalism are necessary. As in^'^'^'O'^^. we consider a 
chain of anyons arranged on a topologically non-trivial 
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cycle wrapping around a torus. Wc denote the space 
of n anyons of types a — (ai, . . . , a„) arranged on such 
a chain with periodic boundary conditions by V^^^^^^^^. 
This space does not naturally decompose into subspaces 
with specified total charge. An orthonormal basis is given 
by the basis vectors 



coefficient of the emtpy graph, i.e. 



|a, b)r 



1 



ffll (22 



-1 In 



bn bi 62 bn-2 bn-l bn 



(14) 



where 6 G fJ" is such that each vertex satisfies the fusion 
rules. The basis specified by Eq. (f4) corresponds to a 
decomposition of the Hilbcrt space as 



^ periodic 



6=(()i,...,6jN 



aibi 



The representation of the bra {a, 6|„ of the vector (14) is 
again obtained by flipping the diagram and reversing the 
arrows, i.e.. 



{a,b\r 



1 bn bi 62 &,i-2 bn-1 bn 

^-p|-....-pp- 

fll "2 In-l In 



(15) 



Local operators acting on a subset of the n anyons are 
represented as before by shaded boxes, but global oper- 
ators can not be represented as superpositions of graphs 
embedded in such a planar surface. Instead, it is con- 
venient to embed the chain and associated anyonic di- 
agrams along one of the fundamental nontrivial cycles 
of the torus. (This is already implicit in Eq. (14).) A 
global operator O : 05ga™ ^pliodk ^ ®aef2" ^pliodk 
mapping between periodic chains of (possibly different) 
lengths TO and n is then represented by a shaded strip 
parallel to the chain on the torus, with m ingoing and 
n outgoing edges, i.e., 

i i I I 

o 



U + 

Application of such an operator is again equivalent to 
attaching it to a diagram. Finally, (partial) traces are 
computed simply by connecting up edges. For example, 
for an operator O acting on 0- Vp^j^j.jQjjij,, the partial trace 
over the n-th anyon is 



tr„ 



O 

i 4 I \ 



TT 



o 



T 



Here we used isotopy on the torus to get a convenient 
expression on the rhs. The (complete) trace is the result 
of connecting up all strands, and then computing the 



tr 



O 



o 



(16) 



Here we denote the coefficient of the empty graph in a 
formal superposition X by [X]vac- That is, it is obtained 
by writing X as a superposition of states with flux a, i.e., 
with a line labeled a going around the torus, and taking 
the coefficient for a = 1. In Appendix A, we show that 
the (partial) trace defined in this fashion is equivalent to 
the orthonormality of the set of vectors (14). 



Anyonic Hamiltonians: long-range effective 
theories 



A remarkable feature of the state space of anyons is its 
topological degeneracy: the Hamiltonian of the quantum 
medium assigns equal energy to each state. Furthermore, 
this degeneracy is stable under local perturbations, a fea- 
ture which makes anyons particularly suited for encoding 
and processing quantum information. These properties 
hold up to exponentially small corrections in the inter- 
anyon distances. 

If the inter-anyon separation falls below a certain 
length scale, the microscopic details of the system be- 
come relevant and the topological degeneracy is generally 
lifted. Such a degeneracy lifting has been examined in 
various quantum media^^^^^. In the system- independent 
anyonic formalism, Bonderson^^ has shown that a gen- 
eral interaction between two anyons can be interpreted 
as tunneling of topological charge, and that generic tun- 
nelling fully lifts the topological degeneracy. 

Hamiltonian terms responsible for such tunneling and 
more generally arbitrary multi-anyon interactions take 
the form of Hermitian operators as in (10) with m = n. 
The exact form of the effective Hamiltonian governing 
the energy splitting depends on the geometric arrange- 
ment of the anyons. In a regular lattice, nearest neigh- 
bor (to = 2) and next-to-nearest neighbor (to = 3) in- 
teractions are most relevant physically as the interaction 
strength decays exponentially with distance. Lattice-like 
arrangments of anyons arise when certain spatial distri- 
butions are energetically favored, e.g., by inserting de- 
fects into the quantum medium which couple to addi- 
tional quantum numbers (such as electric charge) of the 
anyons. 

Paradigmatic models of such effective Hamiltonians 
have been considered extensively in the literature. They 
can be thought of as describing the dynamics of the inter- 
nal degrees of freedom of anyons pinned to fixed sites. We 
refer to^° for an introductory discussion of such models. 
We discuss explicit examples for Fibonacci anyon chains 
in Section V. 
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IV. ANYONIC ENTANGLEMENT 
RENORMALIZATION 

A. The setting 

We have described the origin of anyonie Haniiltonians 
as long-range effective descriptions of quantum media in 
Section IIIC. We now turn to the problem of defining a 
variational ansatz for such systems. For concreteness, 
we restrict our attention to one-dimensional chains of 
anyons, arranged in a periodic fashion along a topolog- 
ically non- trivial cycle on the torus. We discuss more 
general 2-dimensional arrangements in Appendix B. 

One may consider different spatial distributions of 
anyons on the chain. For example, in a setting with sev- 
eral non-trivial anyon types, one may be interested in 
the effective behavior of a staggered, i.e., alternating ar- 
rangement of anyons. While our ansatz could in principle 
be adapted to such cases, here we consider the simplest 
non-trivial setting. We assume that a subset riofi C 
of anyons is allowed in each site. The Hilbert space of a 
(periodic) chain of length n on the torus is then given by 



n 



chain. 1 



eya 
•^pori 



criodic 



(17) 



a (EH", 



tf 



The most commonly considered case (e.g., the golden 
chain^^) is when fieff = {a} consists of a single anyon a, 
that is, each site is occupied by a particle of type a. Our 
formulation is slightly more general, as it allows to con- 
sider Hamiltonians which can create and destroy parti- 
cles on sites of the chain (resp. change particle types) 
when 57ofi = This is important when the quantum 
medium assigns nearly degenerate energies to different 
distributions of anyons. Further intermediate cases could 
be considered. 



B. The ansatz 

Consider a MERA-structure for a periodic chain 
with N sites as in Fig. 1. We associate to this structure 
a family of states in the anyon Hilbert space 'Hcha.in,N 
(cf. (17)) as follows: 

1. To the strip at the top, associate the bra {(p\ of a 
normalized state 



(18) 



2. To every intermediate box with m ingoing and 
n outgoing strands (m > n), associate the ad- 
joint W^^ of an isometric charge-conserving map 



W e End V^, V"' 



(20) 



that is, 



■^.....^ {n strands) 



I-)- 



■f ■f ("i strands) 



3. Regard the state \ip) (Eq. (19)) and the family of 
maps {W} (Eq. (20)) as variational parameters 
specifying a state 1^-) = |^'|^),{w}) & "Hchain.N of 
the chain. The state is determined by the follow- 
ing recipe: in the MERA-structure of Fig. 1, re- 
place every box by the superposition of trivalent 
labeled graphs representing the associated object 
(i.e., Eq. (19) and Eq. (20)). The resuh is a su- 
perposition of labeled graphs, each with N ingoing 
edges. This superposition represents (vPl. 



More explicitly, the maps (20) are of the form 
W = 0W(c) , 



(21) 



that is. 



with W(c) e End(e^^^,,„^ K^ea-'e03, Vc^') satisfying 

W(c)tW(c) = y. . (22) 

Eqs. (21) and (22) severly constrain the set of allowed 
maps W for certain (m, n) and rJeff- 

This ansatz is clearly motivated by entanglement 
renormalization for qudits. In fact, it has the same op- 
erational interpretation: the MERA-structure of Fig. 1, 
after replacing each box by the superposition of graphs 
specified by {W}, is a procedure for successively map- 
ping the chain to a coarse-grained chain by local gates 
and isometrics. Indeed, due to charge conservation, boxes 
with the same number of in- and outgoing edges corre- 
spond to local unitaries on the anyons, while boxes with 
fewer outputs than inputs correspond to local isometrics 
reducing the number of degrees of freedom. Such reduc- 
tions preserve the total charge in their support. Impor- 
tantly, since the range of each operator is contained 
ill ®aGn™ states supported on 'Hchain,Af are mapped 
to coarse-grained chains of the same type at each level in 
Fig 1, i.e., with particles from the subset fioff C 51 on all 
sites. 



where n is the number of ingoing edges. 



(19) 

C. Efficient evaluation of physical quantites 

Having introduced a set of variational ansatz states 
parametrized by (|iy9),{W}), we argue that quantities of 
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physical interest such as expectation values of local ob- 
servables and correlation functions can be efficiently com- 
puted from these parameters. Here we use the formal 
equivalence of the manipulation rules of anyonic states 
and operators with usual tensor contractions. 
The anyonic analog of Eq. ( 1 ) is 



ik Ji ik iL 



^ ^ (n strands) 



(rn strands) ' 



(23) 



for any operator W e End (^e-^^„^ y^ e^-,^,,„ 

with the required isometry property (where m > n). 
Each of the small dark boxes represents the projec- 
tion lo^jj onto the subset of allowed anyon labels. Note 
that, if the edges on the Ihs. of Eq. (23) are connected to 
an operator W as specified in the ansatz, these may be 
omitted. This is because both the support and range of 
W are already restricted to the set of allowed anyons J7off 
(cf. (20)). In particular, we formally recover the rule 
Eq. (1) in this case. 

The second important ingredient is formula (16) for 
the trace of an operator on the chain. Eq. (23) and (16) 
immediately imply that expectation values of local op- 
erators can be efficiently evaluated for a state = 
|^|i^),{w}) in essentially the same manner as for MERA- 
states of qudits. The same is true for two-point correla- 
tion functions for certain distances of the points related 
to the MERA structure. 

Consider for example a local observable O acting on 
two sites of the chain. The expectation value of this 
operator, given a density matrix p, is equal to the dia- 
grammatic expression 



tr(Op) 







di. ii 


1 




■ 1 1 " 









(24) 



If p = |^')(5'| is an anyonic MERA-state corresponding 
e.g., to the structure of Fig. 1 (a), we conclude from (24) 
that 

(*ioi*) = m_ , 

where X is the superposition of trivalent labeled graphs 
specified by the diagram in Fig. 2 (b). Using (23), this 
immediately reduces to [X']^^^, where X' is the super- 
position in Fig. 2 (c). This can be efficiently evaluated 
using the superoperators defined in (3) and (4) (again 
interpreted as anyonic expressions). Two-point correla- 
tions functions can be computed analogously. 



D. Computational cost and refinements of tlie 
ansatz 

To count the number of parameters needed to de- 
scribe the anyonic ansatz states, let ^ = maxagn^jj da be 
the maximal quantum dimension of the particles used. 
Since the number of states of the form (5) can be upper 
bounded by 0(^"), a charge-conserving map as in (20) is 
described by fewer than 0(^"+™) parameters. Similarly, 
a state \ip) as in (18) is described by 0(^") parameters. 

As with the MERA for spin chains, the family of 
ansatz states may be enlarged by replacing an edge 
by s > 1 edges; this is analogous to increasing the 
bond dimension. In this case, isometrics are described 
by 0(^«("+™)) parameters as opposed to 0(#("+")) 
in the qudit case. (Note that, by definition, ^ < 
d for any embedding of anyonic states into qudits: 
for example, for the Fibonacci anyons considered be- 
low, ^ w 1.618.) In summary, a general anyonic 
MERA is described by 0{poly{^^^)N log N) parameters, 
and translation-invariant and scale-invariant MERAs by 
O {poly {^'^) log N) and 0{poly{2!'')) parameters, respec- 
tively. 

The remainder of this story is the same as that of the 
MERA for spin chains, and we refer to the extensive lit- 
erature (e.g.,^^) on this subject. For example, the com- 
putational cost of computing local expectation values is 
roughly the same as for qudit chains (with x replaced 
by ^*), and methods used, e.g., for varying over the 
isometrics in Fig. 1 can directly be adapted to the any- 
onic framework. Compared to the qudit chain setting, 
an additional advantage stems from the fact that the 
isometrics are charge-preserving and thus take a block- 
diagonal form. Therefore, matrix multiplication and sin- 
gular value decompositions can be performed on matrices 
whose dimensions are a constant factor smaller than with 
a naive qudit ansatz. Similarly, methods for extracting 
critical exponents from scale-invariant MERA states^^'^^ 
may be applied in the anyonic setting. 



E. Example: Fibonacci anyons 

The Fibonacci theory has one nontrivial particle r with 

quantum dimension dr ~ (j) ^ ^2'^ equal to the golden 
ratio and fusion rule r x t = 1 + r. The _F-matrix is given 

by 



V. 



' 
77 



(25) 



where we use the convention that a solid line represents 
the r-label, while a dotted line represents the trivial la- 
bel 1. Consider a periodic chain of r-anyons. To get a 
corresponding family of ansatz states, we set rigff = {t}. 
Using the convention that a solid line represents an edge 
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with label r and a dotted line represents an edge with 
label 1, the Hilbert spaces "Hchain.n of a periodic chains 
with n e {1, 2} particles are spanned by (cf. (14)) 



U 



chain, 1 



C|r) 



(26) 



^chain,2 = span | | , | | , ^ |^ 

span{|rT),|lT),|Tl)}. 



In general, the space Hchain,n is spanned by vectors cor- 
responding to (periodic) sequences of Is and ts, with the 
fusion constraint forbidding neighboring Is (this defines 
an embedding into a subspace of (C^)*^"). This deter- 
mines the form of the state \ipn) S "Hchain.n correspond- 
ing to the top box in structures as in Fig. 1 (a) and (b), 
respectively. Next, we consider the constraints on the 
maps W = W(„ in (21) with n input- and m output 
anyons, for small {n,m). The standard form of charge 
preserving isometries/unitaries is 



W(l,2) = 



w 



(2,2) 



01/4 



01/2 



w 



(1,3) 



01/2 



-V. 



^ A 



01/2 



w 



(2,3) 




(27) 



(28) 



(29) 



(30) 



where 9 € [0, 27r[ and |ap -I- = 1. In aU these cases, 
the map is completely specified by a phase and/or a qubit 
state. A less trivial case is W(3_3), which is equal to a 

phase times the projection onto the span of plus a 



two- by- two unitary W(r) on the span 



of{V^,Y}. 



For later reference, we also state the most general form 
of an isometry with two input and four output strands: 



W 



(2,4) 




From (26), (27) and (28), we conclude that the family 
of states associated to Fig. 1 (a) is rather uninteresting 
as the variational parameters {\ip), {W}) are merely a set 
of phases. In contrast, the structure in Fig. 1 (b) gives 
rise to a less trivial family of states due to (29). 

We will give additional non-trivial explicit examples in 
Section V. 



F. Distillable states for composite anyon coding 

The MERA ansatz for qudits is motivated by quantum 
circuits. Indeed, a MERA-description of a state provides 
a circuit preparing the state starting from the top-level 
state \ip). This is achieved by realizing isometrics using 
ancillas prepared in pure states; it corresponds to run- 
ning the coarse-graining procedure in reverse. This use 
of the MERA has been proposed for example as a way of 
efficiently preparing topologically ordered states^*^. 

The one-to-one correspondence between preparation 
circuits of a certain form and states described by the 
entanglement renormalization ansatz clearly extends to 
anyons. However, there is an additional relation in the 
anyonic setting which corresponds to running the coarse- 
graining forwards (instead of backwards as in the qudit 
case): a subclass of the anyonic ansatz is in one-to-one 
correspondence with distillation procedures preparing a 
logical state \ip) of composite anyons. The anyonic renor- 
malization ansatz therefore provides an alternative char- 
acterization of "distillable" states in the framework of 
composite anyon coding'^^. 

The goal of composite anyon coding is to prepare a 
suitable state for computation starting from some un- 
known initial state, but without using measurements. 
In the terminology of anyonic entanglement renormal- 
ization, composite anyons are anyons at higher levels in 
the coarse-graining scheme. Initial states which allow to 
prepare a given target state \tp) of composite anyons can 
be characterized as follows. They can be represented by 
an anyonic renormalization ansatz (corresponding to the 
preparation scheme) with the following properties: the 
state at the top is fixed to \ip), and all adjoints of isome- 
trics W^^ are implementable by braiding and fusion. The 
latter condition means that the renormalization scheme 
only consists of unitaries effected by braiding, and coarse- 
graining, that is, any reduction in the number of anyons, 
is achieved by bunching together some particles. Such 
fusion-based coarse-graining is given by a product of the 
isometry 



(31) ^|a,fe;c)(c|=^ 



dadb 



1/4 



a ~ a(c), h = b{c) 



where I a I 



+ h\' + \s\' 



whose adjoint describes the fusion of a pair of particles. 
An anyonic renormalization ansatz with these proper- 
ties can directly be implemented using the operations 
commonly envisioned to be available for manipulating 
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anyons. Thus we can regard such schemes as state prepa- 
ration circuits for topological quantum computers. 

V. APPLICATION TO THE GOLDEN CHAIN 
AND THE MAJUMDAR-GHOSH CHAIN 

A. The model 

In this section, we consider the use of the anyonic 
MERA ansatz in the context of the golden chain^^ and its 
relatives. An introduction to these models can be found 
in^°. Specifically, we consider the Fibonacci-anyonic ana- 
log of the Heisenberg and the Majumdar-Gosh (MG) spin 
chains. The former was introduced in^^ and consists of a 
uniform chain of N Fibonacci anyons with Hamiltonian 
terms favoring one of either possible total charge of two 
neighboring r-particles. Concretely, the Hamiltonian is 
given by 

JV 

j^. ^golden ^ ^ (32) 

i=l 

where each term H2 is a projection onto trivial charge of 
the anyons i and i-|-l, i.e., it has the diagrammatic repre- 
sentation In analogy with the familiar SU{2) spin 

chains, the case J2 > energetically favoring trivial to- 
tal charge is referred to as 'antiferromagnetic' (AFM) 
coupling, whereas J2 < is called 'ferromagnetic' (FM) 
coupling. 

The Majumdar-Gosh (MG) chain^* is a model 
of SU{2) spin-1/2 particles arranged on a chain, with 
three-particle interactions favoring either total spin 3/2 
(called ferromagnetic/FM) or 1/2 (called antiferromag- 
netic/ AFM). Its anyonic analog^^ takes the form 

N 

J3.ijMG^_j^^^. ^ (33) 
i=l 

where is the projection onto trivial charge of three 
neighboring anyons. Using the F-matrix (25), it is 
straightforward to rewrite the terms H2 and in the 
standard basis (14). Corresponding expressions can 
be found in^^. Expressed in the standard embedding 
into (C^)**", this leads to 3- and 4-qubit terms, respec- 
tively. 

The Hamiltonian (32) was studied in detail in^'. Crit- 
icality and a two-dimensional GET description were es- 
tablished numerically. Furthermore, an exact mapping to 
a standard integrable 2D classical lattice model^^ known 
as the RSOS model was given. These studies were ex- 
tended to the one-parameter family of Hamiltonians 

Hg:= cose- Hg°ldcn j^^^^Q. ffMG ^ (34) 

which exhibits a rich phase diagram as discussed exten- 
sively in^^, see Fig. 3. Subsequent work^^^^^ considered 




MG 



FIG. 3: Phase diagram of the model (34) as obtained in®^. 
Phases I and II are gapped, with the exact ground states 
known at the Majumdar-Gosh (MG) point 9 — 3tt/2 (see 
Section VB) and at tanS = 0/2. There are two extended 
criticial phases for which an exact mapping was established 
to the 3-state Potts and the tricritical Ising model at the FM- 
and AFM-golden chain points, respectively'^^. A small sliver 
of an incommensurate phase is found near 9 = I.OTStt next to 
a phase with Z4-symmetry; both these phases are believed to 
be critical. See''* for a detailed discussion. 

generalizations to su{2)k anyons and established a con- 
nection between the gapless modes of these anyon chains 
to edge states of topological liquids^*, which then pro- 
vides some insight into the collective behavior of anyons 
in a two-dimensional setting^®'^^. We refer to these ref- 
erences for further details, and do not attempt to give 
a complete account here. Instead, we restrict ourselves 
to a few example computations illustrating the power of 
the anyonic MERA ansatz. Specifically, we consider the 
model (34). 

B. Exact RG fixed point at the Majumdar-Gosh 
point 

We first consider the FM (i.e., J3 < 0) case of the MG 
chain, the anyonic analog of the Majumdar-Gosh point 
of the spin-1/2 Heisenberg chain for the point 9 = 37r/2 
in (34). This point lies in a gapped phase extending 
from 6 1.387r to ~ 1.5287r, with four-fold degeneracy 
throughout the phase (for chains with even length) 
The ground space at this point is spanned by the states 

I "\2 J...U...U...U...U...U...L (35) 

ivJ LwJ Lwvl LvwwJ LvwJ Lvwl Lvi 

and their translates by one site. Here, a wiggly line de- 
notes a superposition 

_ _ 1 1 

We will now argue that the two-dimensional subspace 
spanned by (35) is exactly described by a scale-invariant 
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anyonic MERA of a simple form. In other words, these 
states arc fixed points under the corresponding renor- 
mahzation group procedure. This provides an encoding 
of a subspace of the ground space, with the property 
that the identity of the encoded state is revealed only at 
the top level in the MERA structure. This is analogous 
to MERA-descriptions of topologically ordered systems 
which arc exact fixed points of a rcnormalization group 
flow: information encoded in a topologically degenerate 
ground space can be recovered at the top level of the 
MERA-hierarchy^'^'^^. However, in the case considered 
here, there is a local order parameter given by the den- 
sity of T-labels on the chain. 

To specify the scale-invariant MERA-ansatz describ- 
ing the ground space (35), consider the MERA-structure 
of Figure 1(a) with refinement parameter s = 2, i.e., 
with every strand replaced by two. Due to this doubling 
and scale-invariance, the corresponding ansatz then is de- 
scribed by an isometry W = W(2.4) of the form (31) and 
a unitary U acting on four r-anyons. We set the unitary 
equal to the identity, and the isometry equal to 



(36) 



which corresponds to the parameters (a, /?, 7, (5, e) = 
(^,-^,0, ^,^) in (31). This completes the specifica- 
tion of the MERA up to translation (which we fix later), 
as the top-level state \if) depends on the actual state con- 
sidered. Note that the MERA-structure of Figure 1(a) 
with unitarics equal to the identity is known as a tree- 
tensor network^^^''^ in the non-anyonic setting. 

Let us argue that the renormalization group scheme 
has the ground states (35) as fixed points. This is most 
easily seen for the state jlrlrlT • • • ) using the diagram- 
matic calculus: applying a layer of (adjoints of) isome- 
tries corresponds to stacking parallel copies of (36) 
on top of the state. We assume that the isometrics arc 
aligned in such as way that this takes the form 



for the state | Irlrlr • • • ) , where we suppressed scalar fac- 
tors. Here we used isotopy invariance and (7). With (14), 
it is easy to verify that the operation (W^)**^/^ also pre- 
serves the norm of the state. This establishes the claimed 
fixed-point property for | Irlrlr • • • ). To verify the claim 
for the second state \tTxTTj.tTx • • ■ ), it suffices to observe 
that this state is essentially equivalent to the former but 
with r-flux, i.e., 



C. Numerical variation over ansatz states 

To assess the suitability of the anyonic entangle- 
ment renormalization ansatz as a numerical method, we 
have implemented an algorithm for numerically miniz- 
iming the energy by varying over the parameters of a 
translation- invariant ansatz. The algorithm is based on 
iterative optimization of the (identical) isometries at each 
level and the top-level state. It is described in detail 
in^"'^ for non-anyonic spin chains. For a fixed isometry, it 
proceeds by computing its environment, that is, the con- 
traction of the MERA-network with the isometry omit- 
ted. The resulting tensor can be interpreted as a linear 
map whose singular value decomposition dictates how the 
isometry is updated. Adapting this to the anyonic set- 
ting is straightforward: here the environment always has 
a block-diagonal form with respect to total charge. Com- 
pared to the algorithm of^^, the only significant difference 
lies in the implementation of the ascending and descend- 
ing superoperators (see e.g., (3) and (4)) used to compute 
the environments. As with all anyonic operations, they 
require applying basis changes into compatible tree-like 
bases (cf. (6)). These basis changes can be prccomputed. 

This randomized optimization algorithm is suscepti- 
ble to local minima, and its convergence depends on the 
choice of initial points. In practice, these issues appear 
to be minor and can be addressed by starting with a 
large number of initial points and postselecting after a 
few iterations. 



Ground state energy and correlation functions 

We have applied the variational algorithm to periodic 
chains of 12 and 16 Fibonacci anyons governed by the 
Hamiltonian (34). These system sizes were chosen to 
allow for comparison with exact diagonalization data and 
to test the suitability of different MERA-structures. 

For the 12- and 16-anyon chains, we use the "ternary" 
and "binary" MERA-structures 



r 




(37) 



nftj Lwwl Iwwl LvwJ LvwJ Lwuvl__Lw 

because of the identity 

This immediately implies that this state is fixed by the 
same rcnormalization group scheme. 



with s = 2 (i.e., every strand is doubled). The former 
consists of a single level of coarse-graining isometries, 
while the latter has two such levels. The variationally ob- 
tained ground state energy is compared to its exact value 
in Fig. 4. As shown, we find good agreement between the 
variationally estimated ground state energy and its exact 
value, over a wide range of values of the parameter 9. 
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7r/4 7r/2 37r/4 tt 57r/4 37r/2 77r/4 




7r/4 7r/2 37r/4 tt 57r/4 37r/2 77r/4 



(a) Ternary MERA-structuro applied to 12 anyons 




7r/4 7r/2 37r/4 tt 57r/4 37r/2 7-k/A 

9 

(b) Binary MERA-structure applied to 16 anyons 



FIG. 4: Ground state energy approximated variationally by 
an anyonic entanglement renormalization ansatz. The in- 
set in each figure shows the relative error Ai5 = (-Bmera — 
£^o)/(£^max — Eo). Fig. (a) is based on a ternary anyonic 
MERA for a system of 12 anyons, while Fig. (b) shows the 
result of a binary anyonic MERA for a chain of 16 anyons. 
At the Majumdar-Gosh point [6 — 3tt/2), the variational 
procedure recovers the exact fixed-point discussed in Sec- 
tion VB for the binary MERA. In (b), the approximation 
is best around around the golden chain-point with AFM- 
couplings 6^ = 0. This is consistent with the fact that the 
ground state has a Z2-sublattice ordering^^ which is compati- 
ble with the coarse-graining structure of a binary MERA. For 
9 = TT, the ground state has a Zs-sublattice ordering for which 
the binary MERA structure is less suited. However, the ap- 
proximation is still better than in intermediate regions where 
the Hamiltonian has both 2-local and 3-local terms. Simi- 
larly, Fig. (a) shows that a ternary structure appears to be 
suitable for capturing the Zs-sublattice ordering at the golden 
chain-FM-point 6 = tt. 



FIG. 5: The first and second (inset) derivative of the ground 
state energy with respect to the parameter 9 approximately 
reveals the location of the phase boundaries (red vertical 
lines). This data was obtained from the binary MERA for 
16 anyons used to produce Fig. 4 (b). For each of the 40 
points 9o G [0, 27r], the ground state energies Eo{9o ± A6) at 
two neighboring points at distance A9 = 10"'' were approx- 
imated using the MERA. The plots show the corresponding 
discrete approximation to the first and second derivative at 
each point. We stress that only limited information can be 
gained from this plot. In particular, it does not reveal the na- 
ture of the phase transition. We refer to the detailed discus- 
sion in*'', where e.g., the CFT descriptions of the transitions 
out of the tricritial Ising phase have been identified. 



Fig. 5 shows that ground state energies computed us- 
ing the binary MERA are sufficient to obtain a rough 
estimate for the location of the phase boundaries. 

To study whether the anyonic MERA correctly repro- 
duces correlations in the ground state, we have com- 
puted the (translation-averaged) two-point correlation 
functions 



C(r) 



1 ^ 

-Y,{{H'^H^+^)-{H^)-{Hi+^-)) (38) 



of the local topological charge density (as measured in 
terms of the local projection onto trivial charge for 
a pair). The result of this computation for the AFM- 
point 6 — and the FM-point 9 ^ tt arc shown in Fig. 6 
and Fig. 7, respectively. They exhibit a remarkably good 
agreement with the exact values. 

In Fig. 8, we show the error when computing nearest- 
neighbor and long-range correlations as a function of the 
Hamiltonian parameter 9. We observe good accuracy in 
regions where the ground state energy is well approxi- 
mated (compare Fig. 4). As expected, the considered 
MERA-structures are less suited for describing ground 
state correlations at intermediate values of 9, where the 
Hamiltonian has both 2-anyon and 3-anyon interactions. 
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C(r) 
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r 

(a)C(r) with the ternary MERA at the AFM-point 



C{r) 




1 2 3 4 5 6 7 8 



r 

(b)C(r) with the binary MERA at the AFM-point 

FIG. 6: Two-point correlation function C(r) (cf. (38)) of the 
local topological charge density at the AFM point. They 
reveal a Z2-sublattice ordering of the ground state wave 
function. The inset shows the absolute error |(5C(r)| — 

iCMBRA(r) - Cexact(r)|. 



We emphasize that the structures (37) with s ~ 2 
are two of the simplest possible leading to non-trivial 
families of ansatz states for Fibonacci anyons. The 
support of the unitarics/isometries (4 anyons) is only 
marginally larger than that of the terms in the Hamil- 
tonian (34). This suggests that the approximation by 
such states may be rather coarse. Nevertheless, the 
ansatz provides reasonable approximations to the ground 
state energies (for all 9), and correlation functions at 
the AFM- and FM-points where the Hamiltonian con- 
sists of nearest-ncighbor-terms. This illustrates that any- 
onic entanglement renormalization successfully exploits 
the constraints imposed by conservation of topological 
charge. Future work may go beyond this proof of princi- 
ple by considering refined families with parameters s > 2. 
This should lead to significant improvements as in the 
non-anyonic setting: Here accurate results for correla- 
tion functions are usually obtained only for high bond 
dimension (e.g., x = 22 for Ising and Potts chains at 
criticality^^). 



C{r) 




(a)C(r) with the ternary MERA at the FM-point 



C{r) 




r 



(b)C(r) with the binary MERA at the FM-point 

FIG. 7: Two-point correlation function C(r) (cf. (38)) of the 
local topological charge density at the FM point. The inset 
shows the absolute error \5C{r)\ = lCMERA(r) — Cexact(r)|. 

Larger systems 

To test the scalability of the method, we have ad- 
ditionally computed the ground state energy of chains 
of length N e {32, 64, 128} using the binary MERA- 
structure (with s = 2) obtained by adding levels to (37). 
Since this is beyond the reach of exact diagonalization, 
we consider the AFM-point, which allows us to compare 
our results to the CFT-predictions of^*". 

Explicitly, we use the fact that the low-lying spectrum 
of a periodic ID critical quantum systems of length TV 
takes the form 

E = eN+ — (hL + hn--) . (39) 

Here e and v arc non-universal constants, c is the cen- 
tral charge of the CFT, and and hji are the confor- 
mal weights of the holomorphic and antiholomorphic part 
of the local field associated with the energy level. The 
latter parameters are defined in terms of a representa- 
tion of the Virasoro algebra and are tabulated for uni- 
tary minimal CFTs. In^*", the CFT corresponding to the 
AFM-point was unambiguously identified as that describ- 
ing the tricritical Ising model at its critical point, with 
central charge c = 7/10. The ground state energy Eq 
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7r/4 7r/2 37r/4 tt 57r/4 37r/2 77r/4 27r 

e 

(a)|<5C(l)| and |<5C(6)| for the ternary MERA 




7r/4 7r/2 37r/4 tt 57r/4 37r/2 77r/4 27r 

(b)|<5C(l)| and |<5C(8)| for the binary MERA 

FIG. 8: The absolute error \5C{r)\ = |CMERA(r) - C,^^ct{r)\ 
of the nearest-neighbor (r = 1) and long-range (r = 6 re- 
spectively r = 8) charge-charge correlation functions. Quali- 
tatively, these deviations agree with the errors in the ground 
state energy (see Fig. 4). Note that at the Majumdar-Gosh 
point {9 = 37r/2), the plotted \5C\ is large even though the 
MERA accurately represents one of the ground states. This 
is because of the degeneracy and the fact that Ccxact is com- 
puted from the completely mixed state on the ground space. 



corresponds to = Hr = 0, whereas the first excited 
energy Ei is determined by /i^ = /i/j = 3/80. 

Using the exact values of Eq and i^i for = 16, wc 
determine the non-universal constants in (39) (approxi- 
mately). The resulting prediction for the ground state 
energy density Eq/N for system sizes < 16 differs by 
only 10~^% from that obtained in the same way using 
the exact spectrum at = 8. This suggest that finite- 
size effects are negligible, and we use the prediction (39) 
obtained in this fashion to study the anyonic MERA for 
larger systems. 

Fig. 9 shows the result of this computation. We find 
that the ground state energy density is well approximated 
by the anyonic MERA ansatz even for larger systems. 

Throughout this paper, we have considered finite-size 
systems. A modification of the optimization algorithm 
for scale-invariant qudit MERA allows to numerically 
estimate data of the associated CFT in the continuum 
limit^^-^^. Adapting this to the anyonic setting, this can 
provide an additional benchmark when comparing to the 
predictions of^^'^^. Indeed, results of this kind were rc- 
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FIG. 9: Ground state energy per site obtained for A'^ = 
16, 32, 64 and 128 anyons at the AFM-point using the anyonic 
MERA (blue crosses). The red line shows the CFT prediction 
for the ground state energy obtained by extrapolating from 
A'^ = 16 (as discussed in the text). The dotted curves indicate 
the (modulus of the) deviation from this prediction when us- 
ing N = 8 instead. The inset illustrates the relative error of 
the MERA with respect to the CFT prediction. 



cently presented in^^ (see note below) 



VI. CONCLUSIONS 

We have proposed a variational ansatz for periodic 
chains of interacting anyons. It is the natural counterpart 
of entanglement renormalization for spin chains. Based 
on the empirical evidence for the numerical accuracy of 
latter, it is reasonable to expect that the ansatz is a pow- 
erful tool for describing critical anyonic systems. Indeed, 
we have obtained numerical evidence for its suitability by 
comparing with exact diagonalization data in the case of 
(a variation of) the golden chain. 

Our ansatz makes optimal use of the anyonic struc- 
ture of the Hilbert space by incorporating conservation 
of topological charge at different scales. We expect this to 
lead to significant computational savings and improved 
accuracy compared to more conventional methods based 
on embedding the anyons into qubits. It may be more 
pronounced for more general anyon models than the Fi- 
bonacci anyon chains numerically studied here. 

Beyond providing an efficient numerical tool, the pro- 
posed ansatz is a starting point for interesting general- 
izations. For example, along the lines of^""^, one may de- 
fine a Hamiltonian renormalization group flow based on 
the ansatz. This flow generalizes the perturbative renor- 
malization prescription analytically considered In^*-'^'^ in 
the context of random couplings. More importantly, the 
current work may serve as a stepping stone for the devel- 
opment of variational methods for 2D systems of inter- 
acting anyons. We give a rough sketch of corresponding 
adaptations in Appendix B. 

The transition from entanglement renormalization for 
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qudits to anyons is remarkably simple on a conceptual 
level: it boils down to the replacement of isometrics by 
topological charge-preserving isometrics, and a reintcr- 
pretation of networks in terms of the isotopy-invariant 
calculus. This suggest that more generally, tensor net- 
work states such as MPS or PEPS may also be adapted 
to the anyonic setting using similar substitutions. This 
should add significantly to the repertoire of variational 
methods for topologically ordered systems. 



Acknowledgments 

The authors thank Simon Trebst and the anonymous 
referees for numerous useful suggestions. They also thank 
Robert Pfeifcr, Guifre Vidal and Julien Vidal for com- 
ments, and Jens Eisert for bringing*^'^® to their atten- 
tion. RK acknowledges support by the Swiss National 
Science Foundation (SNF) under grant PA00P2-126220. 
He also thanks Station Q for their hospitality and Lukas 
Fidkowski, John Preskill and Ben Reichardt for discus- 
sions. EB is supported by DoE under grant no. DE- 
FG03-92-E40701 and NSF under grant no. PHY-0803371. 

Note added: Shortly after posting this work to the arXiv, 
Pfeifer ct al. have also presented work on the anyonic 
MERA'^2. They use the scale- invariant MERA to nu- 
merically extract GET data in the thermodynamic limit, 
studying the Fibonacci- golden chain at the AFM point. 



where 




Since we are interested in the vacuum coefficient of X, 
we can set all q = 1, getting 




(A2) 



But each 0-like graph on the rhs. is proportional to the 
empty graph, with scalar 



= \/ dadbdc . 



(A3) 



This can be verified by applying (7) twice. Inserting (A3) 
into (A2) gives 

[AT] vac = Y\_ ^bj,b'. J daj ■ 



When combined with (Al), this implies the claim. 



Appendix A: The inner product 
1. Anyons on a torus 

In this appendix, we show that the basis states {|a, 6)} 
of ®a Aperiodic i'^^' (1^)) arc iudccd orthonormal with re- 
spect to the Hilbert-Schmidt inner product defined by 
the trace (16). By definition, we have 



\{<bMb)f = U-j^ 
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Inserting the projection I{(ei,e2)} onto a pair of 
anyons (ei, 62) decomposed into total charge c as 



^{(61.62)} 



ei 62. 



ei 




into the horizontal lines leads to 



(Al) 



2. Anyons on a disc 

For the spaces ®ggs^„ ^.gj^ of n anyons on a disc, the 
(partial) trace can be defined in a similar diagrammatic 
manner as for the space 0- Vp° ^jq^jj^. of a periodic chain. 
This leads to a slight modification when defining the any- 
onic MERA for anyons on a disc. For completeness, we 
include a short description here. 

A partial trace corresponds to connecting up in- and 
outgoing strands of an operator, while simultaneously in- 
serting the operator 




J2dc'\c){c\ 



(A4) 



into each line. (Here we use the convention that either 
input- or output of an operator which acts diagonally on 
anyon labels may be represented by an undirected edge.) 
For example, the partial trace over the n-th anyon for an 
operator U acting on n anyons takes the form 



tr„ 



O 



O 



i i i 
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(a) (b) 

FIG. 10: Two-chain ladder arrangement of Fibonacci anyons 
as discussed in*"', with nearest-neighbor interactions. We con- 
sider the ordering indicated in (b). 



The (complete) trace is obtained by connecting up all 
the strands in this way and computing the vacuum coef- 
ficient, i.e., 



tr 



o 



I + + 




(A5) 



The diagram on the rhs. in (A5) is in fact proportional 
to the vacuum graph. Using (7) repeatedly, it is straight- 
forward to check that the states (5) are orthonormal. 



Appendix B: Braiding and more general 
arrangements of anyons 

In the main text, our focus is on one-dimensional 
chains of anyons. Here we briefly comment on the gen- 
eralization of the anyonic entanglement renormalization 
ansatz to systems of anyons arranged in a more general 
way (e.g., a regular lattice). Such systems have been 
studied before (see^^), and the modifications necessary 
to define corresponding anyonic Hamiltonians are nicely 
explained in*°. In fact, these modifications directly carry 
over to anyonic entanglement renormalization when a 
linear ordering of the anyons at every level is chosen. 
An analogous situation arises when considering fermionic 
tensor networks, and corresponding techniques''''^® can 
thus be extended to anyons. 

For concreteness, we first consider the example of the 
two- leg ladder model discussed in"**^, closely following 
that presentation. This is a system of two chains of 
m anyons each placed next to each other as shown in 
Fig. 10(a) along a non-trivial cycle on the torus. The 
figure also indicates two-anyon nearest neighbor interac- 
tions, which are e.g., projections onto trivial charge as in 
the golden chain. 

The Hilbert space of this system is V^" riodic "^i^l^ ^ = 
(t, . . . , t) , as this space depends only on the topology of 

2m 

the surface (an n = 2m-punctured torus) and the bound- 
ary labels (all equal to t). However, the basis (14) is now 
ambiguous and we must choose a linear ordering of the 



anyons. A convenient choice of an ordering is shown in 
Fig. 10(b). 

Given this ordering, some of the nearest-neighbor 
terms in the Hamiltonian now act on non-neighboring 
pairs (i, i + 2) of strands in the diagrammatic represen- 
tation. To make sense of such terms, it is necessary to 
introduce an additional basis change corresponding to a 
transposition (i, i + i-^ (« + 1, *)• Physically, such basis 
changes correspond to movements of the anyons, and this 
is where the (non-abelian) braid group statistics appears: 
pairs of neighboring anyons may be exchanged in either 
clockwise or anticlockwise fashion. Thus there are two 
inequivalent ways of transposing neighboring pairs. It is 
natural to define the Hamiltonian as the result of trans- 
posing, applying the charge projection, and undoing the 
transposition, averaged over either version of exchanging. 

In the diagrammatic formalism, clockwise- and an- 
ticlockwise exchanges of neighboring anyons are repre- 
sented by over- and undercrossings, respectively. The fol- 
lowing computational rules involving the i?-matrix of the 
modular category are then added to the isotopy-invariant 
calculus: 



b a 



b a 



b a 



b a 



With (6)- (8), these rules imply that over- and undercross- 
ings may be resolved according to 




dadh 





(Bl) 



where the sums are restricted to fusion-compatible la- 
bels c. Eq. (Bl) gives the matrix elements of the linear 
operators Bi,B2 : 0^ V^'' — >■ V^°- corresponding to 
the two types of exchanges. They are unitary if a = 5. 

Given this definition, we can write down the interac- 
tion corresponding to preferred trivial charge of a 'hor- 
izontal' pair (2j, 2j + 2) of neighboring anyons on the 
bottom chain in Fig. 10(a) as follows: 




E 



These are 3-anyon operators acting on the anyons in- 
dexed (2j^ 2j -I- 1, 2j -I- 2). Analogous expressions apply 
to pairs on the upper chain. 

The two-chain ladder is one of the simplest examples 
where transpositions need to be used to define physically 
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interesting Hamiltonians. In a more generally arrang- 
mcnt of anyons, such as a regular 2Z3-lattice, the same 
procedure applies. Depending the chosen standard or- 
dering, however, more transpositions may be required to 
apply a local operator to a subset of anyons. For exam- 
ple, a projection onto trivial charge of two anyons may 
take the form 




where we sum over particles for undirected edges. 

The fact that (B2) is a multi-anyon non-local operator 
now appears to be an obstacle to the use of the anyonic 
entanglement renormalization ansatz. But the special 
structure of (B2) allows to efficiently evaluate its expec- 



tation value for suitable MERA structures in spite of this. 
This is because crossings can be pushed past isometrics 
to higher levels using the fact that 




(and similarly for undercrossings) , and be absorbed into 
the isometries/unitaries by applying the linear maps (Bl) 
to their inputs and outputs, respectively. This eventu- 
ally reduces the evaluation to the expectation value of a 
local operator for a related anyonic MERA-state. The 
efficiency of this procedure depends on the number of 
crossings that need to be resolved, and therefore on a 
judicious choice of orderings. This feature is identical to 
fermionic tensor networks/MERAs, where it is necessary 
to keep track of the number of swaps. 
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